Experimental generation of pseudo bound entanglement 
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We use Nuclear Magnetic Resonance (NMR) to experimentally generate a bound entangled (more 
precisely: pseudo bound entangled) state, i.e. a quantum state which is non-distillable but neverthe- 
less entangled. Our quantum system consists ol three qubits. We characterize the produced state 
via state tomography to show that the created state has a positive partial transposition with respect 
to any bipartite splitting, and we use a witness operator to prove its entanglement. 
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Quantum entanglement has been a source for fascina- 
tion and surprise for more than 70 years [l| . The discov- 
ery that entanglement can serve as a resource for infor- 
mation processing has triggered the development of the 
broad field of quantum information science [i,l3|. Design- 
ing and controlling entangled quantum states has thus 
become a major experimental challenge. Several types of 
entangled states have already been experimentally gen- 
erated and detected in the recent years, e.g. bipartite 
entanglement 

m, tripartite entankementM, , 
well as multipartite entanglement pH [T^. [TslTl^. 
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A particularly interesting_class of entangled states are 
"bound entangled" states [l5j , which carry quantum cor- 
relations of an especially elusive and fragile type. The 
name "bound entanglement" implies an analogy to bound 
energy in thermodynamics, which cannot be freed to per- 
form work. Bound entanglement can be described as 
inherent quantum correlations, without the typical in- 
formation processing potential of "free" entanglement. 
Namely, some correlations of quantum nature are ini- 
tially established during the generation of a bound entan- 
gled state, but nevertheless it is not possible to extract 
a maximally entangled state for any number of copies 
of the state (i.e., entanglement distillation is not possi- 
ble for bound entangled states). However, it has been 
shown that under certain circumstances bound entangle- 
ment can nevertheless be useful to establish a secret key 
p^ . Bound entanglement is expected to exist in Na- 
ture - for example it was shown theoretically in [l^ that 
thermal states of several spin models carry bound entan- 
glement. Bound entangled states are always mixed, and 
due to the small region that they occupy in the space 
of all states, they are vulnerable to decoherence. There- 



fore, it is challenging to create and detect them in the 
laboratory. 

Here, we report on the first experiment in which pseudo 
bound entanglement is created and observed [1^. We 
generate a bound entangled state of a class initially sug- 
gested by Acin et al. [l^ with Nuclear Magnetic Reso- 
nance (NMR), and show their undistillability via state 
tomography, by proving that the partial transposes are 
positive. A suitable entanglement witness is implemented 
to detect the entanglement [1^. As usual in room tem- 
perature NMR we work with states close to the (nor- 
malised) identity, i.e. we consider a state of the form 



l-p 



1+PPbe, 



(1) 



where pbe is the bound entangled density matrix. We 
therefore call the whole state p pseudo bound entangled. 
This is in analogy with the NMR GHZ state [U. NMR 
technology has been widely used to demonstrate quan- 
tum information primitives [H, [13, ISE, HB] , and here 
it proves again to be a versatile tool. 

We consider the family of thrce-qubit states defined in 



PBE = 7V-i( 2|GHZ)(GHZ| 



ai|00l)(00l| +a2|010)(010| 



— |01l)(01l|+ (2) 
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where |GHZ) = ^ (|000) + |111)), and the normalisa- 
tion is = (2-1- Y^^i=i{o-i + The free parameters 
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017^2,03 are real positive numbers that obey the condi- 
tion 010203 ^ 1. This familiy of states is "almost" diag- 
onal, with the only off-diagonal matrix elements coming 
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from the GHZ contribution, see Fig. [2] for a visual repre- 
sentation. The states have rank 7 and are thus far from 
being pure. The states in equation ([2]) have a positive 
partial transposition (PPT) with respect to any bipar- 
tite splitting, but it was shown that they are nevertheless 
entangled [l9l |. Therefore, they are bound entangled. 

In j20| a family of witness operators which detects the 
bound entangled states in equation ^ was presented. 
The family of witnesses has the form 



W = W -el 



(3) 



with 



W = \GHZ-){GHZ 
1 
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(4) 



EyT^(|ooo)(iii| + |iii)(ooo|). 



where |GHZ-) = ^ (|000) - |111)). Note that we are 
not using any specific normalisation of the witness. We 
choose the set of free parameters for the state and wit- 
ness in such a way that the witness detects entangle- 
ment with the highest fraction of white noise. The op- 
timal parameters calculated via convex optimization are 
ai = a2 = az = a = 0.3460 and e = 0.1069 [13]. 

A perfect experimental realisation of the bound en- 
tangled state pbe would lead to the expectation value 
tr(W^PBE) — — £• This is obvious by construction, as W 
is decomposable (i.e. it can be written as the sum of 
a positive operator and the partial transpose of another 
positive operator), and therefore its expectation value 
for a PPT bound entangled state cannot be smaller than 
zero. In fact, ir{WpB¥) = 0. The expectation value for 
the identity part gives — e. The smallest possible expec- 
tation value of the witness W is achieved by the GHZ 
state \GHZ) = -ij (|000) + |111)), because this is the 
eigenvector corresponding to the negative eigenvalue of 
the witness (A = -jr^ - e « -1.03). Therefore, the 



minimal expectation value for a bound entangled state is 
only about 1/10 of the minimal expectation value for a 
free entangled state. This is one of the reasons why the 
detection of bound entanglement is experimentally much 
more challenging than the detection of free entanglement. 

In current liquid NMR experiments a huge ensemble 
of about ~ 10^^ copies of the quantum processors are 
used. The ensemble in thermal equilibrium follows the 
Boltzmann distribution dominated by the spin interac- 
tion with the external magnetic field. The equilibrium 
density operator is well approximated by 




(5) 
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FIG. 1: Schematic carbon NMR spectrum of the CHF group 
of ethyl 2-fluoroacetoacetate and the stick representation of 
the molecule. The J-coupling between the spins as well as 
the H,F spin states in the computational basis are presented. 
We use the spin state representation of the form |C,H,F), 
i.e. C (H,F) corresponds to the first (second, third) qubit, 
respectively. The order of the numbering of the spin states in 
the spectrum is due to the signs of the J-coupling constants. 



where is a spin and system specific constant given by 
Ki — with 7i being the gyromagnetic ratio. Izi 

denotes the z-coniponent of the angular momentum op- 
erator of the i-th nuclear spin and d is the dimension 
of the total system. In usual liquid NMR experiments 
with a magnetic field of ~ 12 T and protons at ambient 
temperatures (T= 290 K) we have kh « 8.4 • 10~^. 

The state of interest in our case is given by eq. H]), 
where we want to generate a specific pseudo bound en- 
tangled state yO, i.e. the deviation from the normalized 
identity, with the probability p available in our system. 
Our 3-qubit system can reach under the above conditions 
a fraction of p for the bound entangled state on the or- 
der of K. To characterize the properties of pbe wc will 
do state tomography to reconstruct the experimentally 
generated state Pcxp- Then we will prove entanglement 
of the part Pcxp.BE- For this purpose we adopt the usual 
witness formalism to the pseudo state case, i.e. we shift 
the witness in eq. ^ by the contribution which comes 
from the identity, namely we introduce the "pseudo wit- 
ness" 



Wnmr = 



W 



P 



(6) 



such that tr(WNMRp) = tr(WpBE)- 

We use the heteronuclear 3-qubit system ethyl 2- 
fluoroacetoacetate [13], consisting of a hydrogen, a car- 
bon and a fiuorine spin (see fig. [T]) with gyromag- 
netic ratios of -/h ~ 26.75 • lO'^ T-^s-'^^jc ~ 6.73 • 
10^T-is-\7i. « 25.18- lO^T-^s-i. 

Using such a heteronuclear system has several advan- 
tages, first we can apply fast local (spin specific) unitary 
operations, because of the large Zeeman energy difference 
between the spins, and second we can focus on the zero 
order coupling Hamiltonian in the rotating frame repre- 
sentation, H = Ji2lzilz2 + Jishihs + J2slz2lzz- Here, 
Ji2 = 161.3 Hz, Ji3 = -190.2 Hz and J23 = 47.0 Hz are 
the j-coupling constants between the spins. 
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We structure the experiment into three parts: 

(i) Generation of a suitable basis state population. 

(ii) Unitary transformations (quantum gates) are ap- 
plied to the diagonal state to generate pbe- 

(iii) The generated state is characterized by state to- 
mography and witnesses. 

In the first part of the experiment we generate a spe- 
cific diagonal state, i.e. 



1 p 
o o 

dcIzlIzS + ddlz2lz3 + delzllzlizz) 



(7) 



with the coefficients dn = ^^'^ , ~ —0.78, dh — 

" a(3a+2)+3 ' " 



2(a+l)^ 
' a(3a+2)+3 
48 



—0.21, dc = 2da, dd — 2di, and = 

a(3a-H2)+3 " ~ 3.85. We scalcd the traceless spin op- 
erators in such a way that the scaling factor p of this 
deviation density operator corresponds to the probabil- 
ity of our pseudo bound entangled state of eq. ((1]). 

Five experiments are performed to generate different 
basis state populations. The five diagonal states have 
the form 



Pi 

P2 
P3 

Pi 

P5 



1 Kff 

8 8 

1 KH 

8 8 

1 KH 

8 8 

1 KH 

8 8 



(3.77/,i/,2/.3) 

(-2/.1/.2) 

(-1.88/,i/,3) 

(-2/.2/.3) 



(8) 



I+'y- ("^-1 - 0-2742 - 0.27 7,3) . 



The manipulation of the equilibrium states to generate 
each of these states is done by standard hard spin se- 
lective pulses, j-coupling evolution [H and at the end 
^-gradients to destroy the unwanted a:, y magnetization 
terms, i.e. we apply spatial averaging [28j . 

The states of these experiments are added with ap- 
propriate probabilities to achieve our dia gon al density 
operator pd in eq. (O (temporal averaging j29l. [soj ) . i.e. 



Pd 



(9) 



qi w 0.36, 92 = 95 ~ 0.27, 93 ~ 0.29, 94 « 0.08. 

It follows ff^. 

On the initial state in eq. ([9]) we apply a set of uni- 
tary transformations to produce the pseudo bound en- 
tangled state pbe- In detail we first apply an effectively 
line selective — -l-rotation around the y-axis in the space 
1 000) , 1 100) without affecting the remaining basis states. 
In the following set of hard rf-pulses and j-coupling evo- 
lutions two gates similar to controUed-NOT gates (with 



qubit 1 as control, qubit 2 as target and qubit 3 as target, 
respectively) are applied |21j] . The total unitary transfor- 
mation of these quantum gates is 

U =^ (|000)(000| + |000)(100| - |lll)(000| + |111)(100| 
v2 

-hi|00l)(00l| -i|010)(010| -f |01l)(01l| 
+ |l00)(lll| -h«|lOl)(llO| -i|llO)(lOl|. 

(10) 

The generated state is characterized by state tomog- 
raphy in the usual NMR setting [29|, i.e. before the 
detection the spin selective rotations Y1E2E3, E1E2Y3, 
E1E2X3, Y1Y2E3, E1X2X3, Y1Y2Y3 and X1X2X3 are 
performed. Here X^ (Y^) denotes a 7r/2 rotation of the 
i-th spin about the x- (y-) axis and Ej is the identity op- 
eration. In the experiments only the carbon spin (see fig. 
[U is detected. To get the information of the H (F) spins 
an appropriate swap gate to the carbon and H (F) spin 
is performed before detection. The experimental data 
gained by this procedure leads to an overdetermined set 
of linear equations for reconstructing the density matrix. 
We reconstruct the density operator from the experimen- 
tal data by a least squares fit (32 ]. 

The experimentally generated state is shown in Fig. 
[2l The excellent agreement between ideal and generated 
state reflects the high controllability as well as low de- 
coherence of our NMR system. The main deviations are 
due to partial relaxation and pulse imperfections during 
the experiment [4l| . 

The fraction p of the pseudo bound entangled state in 

Pcxp — ^^g^'' 1 +PPBE is in our case about p « 10~^. 

How close this state is in comparison to the ideal one 
can be expressed in terms of the overlap or the distance 
between real and ideal state. We find a high Uhlmann 
fidelity [s^l for the pseudo bound entangled state, namely 



Fu = tr^y^ y^BE,thoo/3BE,cxp-\/PBE,thoo^ = 0.98, (11) 

and a small trace distance (33 |. 

dt = = 0.09, with X = PBE,thoo - PBE.cxp- 

(12) 

The standard deviation of the density operator elements 
are obtained from the least square method. They are 
used to calculate the standard deviation of the witness 
expectation value by error propagation. This deviation 
accounts for the inconsistencies in the overhead of the 
experimental data, which are due to imperfections in the 
tomography procedure. The reconstructed state /Obe,cxp 
has a positive partial transpose with respect to any bi- 
partite splitting. The witness expectation value of the 
experimental state is 

tr(I¥pBE,oxp) = -0.029 ± 0.010, (13) 

i.e. the state is bound entangled. The ideal expectation 
value would have been —0.11. Remember that the wit- 
ness expectation value can vary from —1.0 to -f 1.8, i.e. 



4 




the possible range —0.11 < (W) < for bound entan- 
glement is rather small - thus, already small deviations 
from the ideal situation can have large effects. 

In conclusion, we have generated and detected pseudo 
bound entanglement with an NMR experiment. We 
reached a high fidelity, F = 0.98, of the created state 
with respect to the ideal one. The partial transpositions 
of the state were shown to be positive with respect to any 
bipartite splitting. The pseudo witness that we used had 
a negative expectation value, thus proving that the state 
was indeed entangled. Here we calculated the expecta- 
tion value of the witness operator using the reconstructed 



state. This is more appropriate than to detect the wit- 
ness directly [sPl, because some experimental errors are 
"averaged" during the least square method used in the 
state tomography procedure. - Thus, we have shown that 
bound entanglement can exist in the laboratory. 

After completion of this manus crip t we learned about 
related work by Amselem et al. [4(j |. 

We thank Otfried Giihne and Matthias Kleinmann for 
enlightening discussions. We acknowledge financial sup- 
port from the EU Integrated projects SCALA and SEC- 
OQC, as well as from Deutsche Forschungsgemeinschaft 
(DFG). 
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